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The diquark condensate susceptibility in neutral color superconductor at moderate baryon density 
is calculated in the frame of two flavor Nambu-Jona-Lasinio model. When color chemical potential 
is introduced to keep charge neutrality, the diquark condensate susceptibility is negative in the 
directions without diquark condensate in color space, which may be regarded as a signal of the 
instability of the conventional ground state with only diquark condensate in the color 3 direction. 
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It is generally accepted that, the cold and dense 
quark matter favors the formation of diquark conden- 
sate and is in the superconducting phase 1], which may 
be realized in compact stars and, in very optimistic 
cases, even in heavy-ion collisions. To have a stable 
and macroscopic color superconductor, one should take 
into account the electric and color charge neutrality 
condition 0, Q which leads to a new phase, the gap- 
less color superconductivity Q or the breached pairing 
phased . In this new phase, the most probable tempera- 
ture is finite but not zero[6|, and the Meissner screening 
masses squared can be negative 0- In two flavor case, 
the color neutrality can be satisfied by introducing a color 
chemical potential /ig in the four-fermion interaction the- 
ory at moderate baryon density |3j, or by a dynamic gen- 
eration of a condensation of gluon field A® in the frame 
of perturbative QCD at extremely high baryon density|| 
where the back ground gluon field (A®) plays the role of 
the color chemical potential /Lis- 

In conventional approach to investigating two flavor 
color superconductivity, one chooses the first two col- 
ors participating in the diquark condensate and let the 
third one be free:9j. When the Lagrangian in the study 
is color SU (3) symmetric, the ground states with differ- 
ent color breaking directions can be connected by some 
SU(3) transformations in color space and hence consid- 
ering only one direction is enough to describe the system. 
However, when the color-charge neutrality is taken into 
account, one has to introduce nonzero color chemical po- 
tentials at moderate baryon density in the Nambu-Jona- 
Lasinio(NJL) type model, and the global color SU(3) 
symmetry of the NJL model is then explicitly broken. 
In this case the state with only two colors participating 
in the Cooper pairing is probably not the ground state, 
and one has to calculate thermodynamical potentials for 
different states and chooses the one with the minimum 
thermodynamical potential as the ground state. In this 
paper, we will calculate the diquark condensate suscep- 
tibility in the other two directions without diquark con- 
densate in color space, from which we can judge if the 
conventional ground state is stable when the color charge 
neutrality in considered. 

We work in the widely used NJL model applied 
to quarks [lfj. The model has been successfully used 
to study chiral symmetry restoration|lOt lll| . isospin 



symmetry spontaneously broken |l2j Il3l. and color 
superconductivity 0, H 0, El El EM3H E3 at mod- 
erate baryon density. In chiral limit, the flavor SU(2) La- 
grangian density including quark-quark interaction sec- 
tor is denned as 

+ G D (^ Qi7 V^^) (^V^V^) , (i) 

where Gs and Gd are, respectively, coupling constants 
in color singlet and anti-triplet channels, ip c = Gvb T and 
ijj c = vb T C are charge-conjugate spinors, C — ij 2 -f° is 
the charge conjugation matrix, the quark field ipi a with 
flavor index i and color index a is a flavor doublet and 
color triplet as well as a four-component Dirac spinor, 
t = (ri, T2, T3) are Pauli matrices in flavor space, and 
and e"^ 7 are, respectively, totally antisymmetric tensors 
in flavor and color spaces. We focus in the following on 
the color symmetry breaking phase with nonzero diquark 
condensates defined as 

A 7 = -2Gr,<# Q i 7 V' e Q ^/5>, 7 = 1,2,3 . (2) 

To simplify the calculation, we assume that the chiral 
symmetry is already restored in this phase. This assump- 
tion is confirmed when the coupling constant Gd in the 
diquark channel is not too large 15J. To ensure color 
and electric neutralities, one should introduce a set of 
color chemical potentials fi a (a — 1,2, ...,8) with respect 
to color charges Q\, Q2, Qs and an electric chemical 
potential /i e with respect to the electric charge (electrons 
are included) . The ground state of the system is deter- 
mined by the gap equations which can be obtained by 
minimizing the thermodynamical potential f2, 

dn/dA 7 = , (3) 

and the charge neutrality condition, 

Q e = -dCl/dfle = 0, 

Q a = -aO/9// o = 0, (o = l,..., 8). (4) 

The above twelve coupled equations © and (0} deter- 
mine self-consistently the physical condensates A 7 and 
the chemical potentials \i e and fii as functions of temper- 
ature T and baryon chemical potential fib- 



2 



We first consider the conventional ground state with 

A 1 = A 2 = 0, A 3 = A^0. (5) 

In this case, Qi, • • • , Q7 vanishes automatically, and only 
Qs 7^ 0. Therefore, we can only introduce the color chem- 
ical potential /Lt 8 with respect to the color charge Qs to 
ensure color neutrality, and the quark chemical potential 
matrix 



fl = diag((l u i, (1 U 2, [1 U 3, Vdl, fJ>d2i fJ-d3) 



(6) 



in color and flavor space can be expressed in terms of the 
baryon chemical potential electrical chemical poten- 
tial fi e , and color chemical potentials us, 

Mui = Mu2 = Mb/3 - 2/i e /3 + ^ 8 /3 , 

Mu3 = Mb/3 - 2/i e /3 - 2^i 8 /3 , 

Mdi = Md2 = Mb/3 + ^ e /3 + /x 8 /3 , 

Md3 = Mb/3 + Me/3 - 2^i 8 /3 . (7) 

Generally, ^ 8 is nonzero, and the color symmetry 
of the NJL Lagrangian is explicitly broken down to 
SU{2) (g) U(l) with generators Ti,T 2 ,T 3 and T 8 . In 
mean field approximation the thermodynamical poten- 
tial f^o of the system can be evaluated as 



q 



A 2 T 



4Gr 



(S) 



where the last term is the contribution from the elec- 
tron gas. In the modified 12-dimensional Nambu-Gorkov 
space with color and flavor indices, defined by the field 
vector 



* = (V'ul, 1Ph> Vto, Ipdl, i>u1i 

$u2, Ipdl, $u3, i>la, 1pd3, ^3) ) 



(9) 



the mean field quark propagator 5 m y is diagonal and can 
be expressed as 

^mf ~ di a g ( S u \^d2 ^d2Ail S(tl,u2 &u2,dl &u3,u3 ^d3,d3 ) 

(10) 

where the diagonal blocks are defined as 



CX _ ( [#L 



H5X 



with the free quark propagators 

[Go] ia = (iUnJO - P • 7 ± M»a7o) 



(11) 



(12) 



Since the quark propagator is diagonal, we can 
analytically take its trace in the color, flavor and Dirac 
spaces and make the Matsubara frequency summation, 
the thermodynamic potential can be expressed as an ex- 
plicit function of A, /i e , ua, fib and T. Minimizing the 
known Qq, we then obtain 3] in the color symmetry spon- 
taneously breaking phase the gap equation 

uft/^y'-W-^-L, (13) 



(2tt) 3 ^ 



E\ 



and the charge neutrality condition 



Qe = 



e=± ""A 

(/(^ 3 ) + /(^ 3 ))] = 0, 

/(^E[^(i-/^)-/^)) 

3 (/(25« ) - f{EL)) - e (2/(^3) - /(£&)) 
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m: 

3tt 2 







(14) 



where the quasi-particle energies are defined as E^ — 
El T <5m, E± = V(|p|±M) 2 + A 2 , £± = | P | ± p,, E± - 
|p| ± Mu3 and E^ 3 — |p| ± M<i3 with the two effective 
chemical potentials M and S/i given by M = Mb/3 — /i e /6 + 
Ms/3 and 8fi = Me/2, and /(a;) = 1/ [e x l T + l) is the 
Fermi-Dirac distribution function. The equations i|13|) 
and l|14jl determine simultaneously the order parameter 
A and chemical potentials /x e and /Lt 8 in the conventional 
ground state of neutral color superconductor. 

Now we come to the question whether the conventional 
ground state defined through (J5J is stable. To answer this 
question we calculate the second order derivation of the 
thermodynamical potential Q with respect to the diquark 
condensates Ai and A 2 , 



2 Q 



d 2 Q 



Ai=A 2 =0 



dA 2 



Ai=A 2 =0 



(15) 



We call the quantity n diquark condensate susceptibility. 

Since the explicit form of Q with finite condensates Ai 
and A 2 is not easy to obtained, we shall use the method 
of perturbation. We take 



Ai = 5i < A , 
A 2 = 6 2 < A , 

and apply the Taylor expansion 



Qs - Qq = (^1 



2 ' S 2 ) 



(16) 



(17) 



Here the perturbed thermodynamical potential Qs can 
be written as 



Qs = 



A 2 + 61 + S 2 
T ^ [ d 3 p 



M< ; 



12tt 2 



Tr In 



2 ^ J (2tt) 3 
with the matrix T 62) defined as 



S-)+T(S 1 ,d2) (18) 
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(19) 



Taking into account the expansion 

00 (—i\ n 

TrlniA- 1 +B) = TrhiA- 1 - V ^ — >—Tr(AB) n (20) 



thus one can choose Ai = A2 — 0, A3 ^ without los- 
ing generality. When explicitly considering /ig 7^ 0, the 
function K(x) is related to the color charge density Q% 
through [I3 



for general matrices A and B and the relation 
Tr[S mf T(5i,S 2 )} =0 , 



(21) 



which ensures the linear term to be zero, we can expand 
Q$ to the quadratic terms in <5i and S 2 , and obtain the 
diquark condensate susceptibility 



1 Ty, 

2G D + 4 ^ 



"" P -Tr(5 m/ r(l,0)) 2 



(2tt)3- 



(22) 



The next task is to calculate n based on the known 
A, fi e and fi% in the conventional ground state. After a 
somewhat complicated but straightforward algebra cal- 
culation and using the gap equation (|13[1 . k can be sim- 
plified as 



k = -2fj, s K(0) 



(23) 



with the function K(x) given by 



1 \ f(K 3 ) + f(E&) - 1 



£75 



J_ M f{E e +) + f{EL)-l 
where the x-dependence is hidden in F 12 defined as 



(24) 



■ Ms T Eq =f E^ 

■ Ms T Eq ± 



(25) 



It is clear that, in the case without considering the 
color chemical potential, namely /is = 0, we have k = 0, 
the conventional ground state is stable under the pertur- 
bation. In fact, when /ig = the thermodynamic po- 
tential depends only on the quantity ^ + A?, + A§, 



K(-Hs) = = 



(26) 



we can then expand K(x) around x = —/is- Since the 
magnitude of ua is only a few MeV in the whole color 
breaking phase[3Lll8l|. which is much less than the order 
parameter A, the baryon chemical potential fib, and the 
momentum cutoff A, we can keep only the linear term in 
the expansion, 



4/i| f d 3 p ^ r(£^) 2 + (E%) 



E\ 



< - 



A " J ( 2 ^) e"=± ^ 

(1 - f{E%) f(EL)) - 2EI (1 - f(E* u3 ) - f{E%)) 
4/xg f d 3 p ^ { {E^-{El u +El d )/2f 



4^8 f d J p r 
A 4 J (2tt) 3 I 

(l-f(E' + )-f(El))] 



E\ 



< 



(27) 



Therefore, we have proven analytically that, under 
the condition |/i 8 | <C A. fib, A, n is negative in the con- 
ventional ground state when the color charge neutral- 
ity is considered. This conclusion can also be proven 
numerically. To do numerical calculation, we choose 
the parameters G s = 5MGeV 2 , A = 0.653GeV and 
Go = 3G5/4 as in Ref.@]. With these parameters the 
minimum baryon chemical potential where the color su- 
perconductivity phase starts is fib I '3 = 330 MeV. The 
baryon chemical potential dependence of fig and k, cal- 
culated through JT3J), lO and is shown in Fig[T] 
in the case with both color and electrical chemical po- 
tentials and in Fig|21 with only color chemical potential 
(namely taking /i e = 0). We see that, k is negative in 
the whole color breaking phase, independent of whether 
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FIG. 1: The color chemical potential us and k defined in (12.'-H 
as functions of baryon chemical potential /it, at temperature 
T = for color superconductor with both constraints Qg = 
and Q e = 0. 




350 400 450 500 550 600 350 400 450 500 550 600 



/i b /3 (MeV) rt,/3 (MeV) 

FIG. 2: The color chemical potential fj.% and ft defined in I23H 
as functions of baryon chemical potential /if, at temperature 
T = for color superconductor with only the constraint Q$ = 
0. 



the electrical charge neutrality is taken into account or 
not. In fact, in both cases with and without electrical 
charge neutrality, fis is really very small, compared with 
A, /if,, A which are all the order of hundreds MeV, the 
assumption used in the analytic calculation is safe. 

In summary, we have investigated the relation between 
the explicit color symmetry breaking and the ground 
state which reflects the spontaneous color symmetry 
breaking. For color superconductor without color neu- 
trality, the Lagrangian of the system is color SU(3) sym- 
metrical, and the spontaneous breaking can happen in 
only one of the three directions in color space, namely 
from SU(3) to SU(2). However, when the color chem- 
ical potential is taken into account, the Lagrangian of 
NJL model loses its color SU(3) symmetry explicitly, and 
therefore the spontaneous color breaking in only one di- 
rection (A3) may be impossible. Through calculating 
the diquark condensate susceptibility n in the other two 
directions, we have proven analytically and numerically 
that k is always negative in the whole color breaking 
phase in the frame of flavor SU(2) NJL model. This may 
be a signal that the conventional ground state in unsta- 
ble. More detail investigations in this aspect will be done 
in our future works. 
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